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1.  INTRODUCTION 


This  report  describes  an  analytical  method  and  its  associated  computer  program  which  can  be 
used  to  predict  helicopter  vibrations  for  undamaged  and  damaged  main  rotor  blades.  The  analysis  is 
generic  and  can  be  used  to  analyze  any  helicopter  with  a  minimal  amount  of  physical  input  data. 

The  model  consists  of  independent  main  rotor  blade  flapping-torsional  degrees  of  freedom  and 
helicopter  longitudinal-vertical  degrees  of  freedom.  Level  flight  trimmed  conditions  are  determined 
first,  after  which  damage  is  imposed  onto  the  main  rotor  blades.  Vibration  levels  are  calculated  in  the 
cockpit  both  for  the  undamaged  and  damaged  cases.  The  model  considers  aerodynamic  forces  and 
moments  using  blade  element  theory  wherein  each  harmonic  is  calculated  individually. 

This  program  calculates  cockpit  vibrations  for  a  bare  helicopter  with  no  vibration  absorber 
devices  present,  and,  in  conjunction  with  input  of  human  tolerance  to  vibration  data,  the  code  can  be 
used  to  predict  aircraft  kill  categories. 

There  are  other  sources  contributing  to  helicopter  vulnerability  due  to  main  rotor  blade  damage 
besides  human  vibration  tolerance,  which  can  also  place  a  helicopter  in  jeopardy.  These  include  loss 
of  main  rotor  blade  structural  integrity,  vibrations  transmitted  to  other  critical  structural  ccnnponents  on 
the  aircraft,  and  lack  of  pilot  control  authority  to  compensate  for  blade  damage.  These  additional 
factors  are  not  included  in  this  analysis. 

The  aerodynamic  terms  are  developied  using  blade  element  theory  in  Section  3  and  are  referenced 
to  Figure  1.  The  basic  theory  is  from  Barnes  (1967)  and  Chopra  (1992).  In  this  analysis,  the 
aerodynamic  forcing  functions  are  separated  into  individual  harmonics  so  as  to  be  able  to  trace  the 
frequency  filtering  effect  from  the  rotating  blade  into  the  fixed  (fuselage)  system  and  to  isolate  and 
identify  the  fundamental  physical  parameter  contributions  to  the  aerodynamic  forcing  functions. 

Section  4  describes  results  of  application  of  the  analysis  to  a  UH-60A  helicopter  in  which  one 
blade  in  the  main  rotor  is  considered  to  be  damaged.  The  damage  consists  of  the  removal  of  various 
amounts  of  the  outboard  end  (tip)  portions  of  one  rotor  blade.  Section  S  discusses  the  results  of  the 
damaged  rotor  blade  analyses,  and  Section  6  draws  conclusions  from  the  results. 
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Horizontal 


Figure  1.  Aerodynamic  airfoil  definitions. 
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The  Appendix  contains  a  description  of  the  code,  the  input  variables  required,  a  description  of  the 
computer  code  output,  and  a  list  of  output  variables. 


2.  NOMENCLATURE 


a  Angle  of  attack 

a  Lift  slope 

P  Flap  angle 

c  Chord 

Cj  Blade  drag  coefficient 

C|  Blade  lift  coefficient 

C„  Blade  torsional  moment  coefficient 

L  Lift 

D  Drag 

e.  Blade  torsional  elastic  axis, 

+  measured  from  PA  toward  leading  edge 
F  Force  parallel  to  rotor  shaft 

Fq  Steady  value  of  F 

Y  Rotor  shaft  built  in  tilt  angle  +  forward 
Q  Rotor  speed 

P  Helicopter  pitch  angle  +  into  the  wind 

PA  Blade  pitch  axis,  +  measured  from  leading  edge 
n  Pi  =  3.14 

'P  Rotor  azimuth  angle 

r  Variable  blade  radius 

R  Blade  radius 

p  Air  density 

t  Time 

0  Pitch  angle 

U  Average  downwash  velocity 

X  Helicopter  longitudinal  displacement 

Z  Helicopter  vertical  displacement 

N  Number  of  blades  per  rotor 

Q  Coefficient  of  in  C^  equation 

K  Harmonic  Number 
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3.  DERIVATION  OF  AERODYNAMIC  TERMS 


The  rotor  blade  lift  and  drag  per  unit  length  is  defined  as* 


it-l  pC,c(Ur^*Vfh 

dr  2 


where 


+  Up^). 

dr  2 


Cj  =  aa,  a  =  6  -  4> , 


(t>  =  TAN —L  , 

Ur 


Ur^  ar+(£LcOS(P*y)  -  —  SIN (P  * y) SIN 'P 
'  dt  dt 


Up  =  U+  ^SINiP-^y)  *  ^COS(P*y)  +  , 

dt  dt  dt 


0  =  E  COSi^  +  e,,  SIN 


P  =  £  Pi,  C05i'P  +  P„SW<'P 


Resolving  dUdr  and  dDIdr  parallel  to  the  rotor  shaft. 


—  =  —  COS^  -££  SIN^ 
dr  dr  dr 


*  Equations  1  to  4  are  standard  blade  element  theory  equations  from  Barnes  (1967)  and  Chopra  (1992).  The  equations  5  to 
42  were  developed  by  the  author. 
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For  high-speed  level  flight,  <J)  =  U pi Uj,  then 


=  Ipac  {up*Up)  [b-HI. 
dr  2^  '  ^  Ur 


Also  for  level  flight,  let  C(?S<j»  =  1  and  SINfy  =  <t),  so  that 


where  terms  of  third  order  are  neglected 


^  =  [e  [u^  *  uf\-Up  Uj\-^pCoC 


UpUr, 


Collecting  terms 


^  =  ^pc  ae  [uj  +  Up]  -  Uj  Up  (a  +  C^) 


Uj  and  Up  are  written  as 


Ur  =  ar  +  \—-^iP-^y)]siN  ^ 
dt  dt 


TI  TJ  dZ  rdp 

Up  =  f/  +  _(/’  +  Y)+_+_^  . 
dt  dt  dt 


The  downwash  velocity  is  calculated  as 


2pn«' 
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For  level  high-speed  flight.  Up  <  <  Uj,  so  Equation  12  can  be  written  as 


—  =  Ipc  aQUp  -  Up  Up  (a  +  C^) 


Define  the  following  parameters: 


AU  ^  ^ 

dt  dt 
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DU  =  2Clr(AU), 
Eu= 

2 


GU=  U  +i£(F+Y)  +—  . 

dt  dt 


Uf  =  (BU)  +  (DU)SIN^  +  (EU)COS2^ 


Up  Up  =  nr(GU)  +  +  (AU){GU)SItPV  +  (AU)r^SIfPV  . 


Performing  a  radial  integration  of  the  terms  containing  r,  evaluating  the  trigonometric  product  terms, 
and  substituting  in  Equation  16,  we  obtain  new  harmonic  terms  that  are  sums  and  differences  of  the 
original  harmonic  products.  This  is  called  harmonic  cross  coupling. 


This  allows 


if  p  I  ^ 

F  =  4pc  aQ}^^{AU)^^  J^e^.COSi'V 
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3  D  ^ 

+  0  *  {auY—  53  e„  sin 

L  3  2  1 


+  ailR^^^El  530,,.  [S/^(/  +  l)'F-S/N(i-l)'P] 


2  • 


+  0„  [C0S(i-l)4'-C0S(/  +  l)'J^] 

2 


-  fl  (>1{/)2£5^  [^^<55(i+2)  y  +  COS(j-2)4'] 


A 

a  {Au9—  53  0.,  [5/yV(/  +  2)‘P  +  S/N(i-2)4'] 

4  i 


(fl  +  Co)£/?2(Gf/) 


(a*Co)R^^Yl  «(-P„S/Ni'F  +  Pi^COSiT) 


-  (a^Cp)(AU){GU)R  S//V'P 

p2^-lR.f 

-  {a-^Cp){AU)!L.aY,  _ —  [COSO  -DT-CGSCj +  1)'P] 

2  imo  2 


(a  +  Cp)(>lt/)^n53l^  [S//V(j  +  l)'P-5/Ar(/-l)4']  . 

2  imo  2 


For  torsion,  the  moment  per  unit  span  is 


dM  1  ^  ^2  I  rr2  ,,2\ 
__  *  _  P  C  [(^T  /  ' 


where 
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(26) 


Cm  *  CMo*ycL*Qa^  . 


r  ^ 


Cm  “  C^fQ  +  Y 


0-. 


f/i 


r  r,  ^2 


0-. 


t// 


and 


=  Ipc^  (c^o  +  Y0(/^  -yUjUp  +  G0^i/r  -  IQ^UjUp] 
dr  2  ' 


After  integrating  radially  along  the  blade  and  separating  out  harmonics  similar  to  the  harmonic 
coupling  that  was  obtained  for  the  force  terms  F  of  Equation  24. 


M  =  Lpc^  Cj^fQ 


>2/? 


LV. 


3  2 


+  aR^(AU)SIN'¥-iAU)^lcOS2'¥ 


+  ^pc^Y 


aR 


aR 


LV 


3  2 


5:0..,COSi'F 


l«0 


^ipc^r 


,2/? 


aR 


£2^  _+(/!(/) 

3  2 


aR'^l^El  ^0,^  [S/1V(j  +  1)'P-S//V(i-1)'F] 


1  2 
+  _pc^Y 


[Ir2(AU1  ^0.^  (C0S(i-l)'F-C05(i  +  l)'F] 

2  f 


V 


D  " 

y;  0,JcOS(i  +  2)'F  +  C(9S(i-2)'F] 
4  ,-o 


(27) 


(28) 

cross 
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^pc^Y 


(A U)^  1  53  e,,  [sin  ( j  +  2)  4'  +  5/iV (I  -  2)  'F] 


4 

l-O 


^  K  ^ 

4pc2Y|«^(G£/)+/?3ii!  53/(-p..,5/A^/4'^P,,C(9S/4') 


1*0 


1pc2y(A{/)(G(/)/?S//V4' 


Ipc2Y(AG)/?2ii  ([C0S(i-1)4'-C05(j  +  1)‘F]) 

2  2  2 


-^pc^y(.AU)R^E'£  i^i,lsIN(i*\)^-SIN(i-\)'¥] 


*jpc^Q 


>2R 


.2R 


£2^_  +  (Af/)‘‘_ 
3  2 


^  0^ 

53  _li  (cos  2i4'  +  l) 
i^o  2 


K 


53  0.,0.,SW  2/4' 


»2/? 


.2/? 


fl^_  +  (A(/r_ 
3  2 


^  0^ 

53  _li  (l  -cos  2/4') 

■  2 


+  ^pc2Gn/?^l^^^  o2  f  1 


a:  / 

53  QJ^  1  [s//V(2/  +  l)4'-S//V(2/-  l)4']+S/W4' 

^■o  V  ^ 


+  Ipc2Qi2/?2ll^  53  0,^^  [cos  (2/-l)4'-COS(2/  +  1)4'] 


1*0  ^ 


+  Ipc2en/?2^£G)  02  r5/;v*p_2  [s//V  <2/  +  1)4'-S/Ar  (2/-1  )4']'' 

2  2  2  y 
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^  'N 

E  ^fc  lCOSi2i*2)H>  *COS(2i-2)'¥]*COS2H‘ 


4  ^  "^12 
^1*0 


ipc^GOlOT^i  E 

^  ^  i^o  ^ 


^  0.,0. 


[s//V(2j +2)4'+S/A/(2/-2)4'] 


ipc^ec/ii/)^*  £ef/ 

^  I  B/l 


C(952'F-lc05(2i  +  2)'P+C0S(2j-2)‘F 
2 


lpc^2Q^i^(GU)  X:  (0.,COSi4'*0..,S/^i'F) 
2  2 


-Ipc22QJ22^  I:  0.,f:  -^[SINii*p)^-SIN(i-p)'¥] 


^  ic  2 


-lpc^2Qn^^  £  0.,f:  £^[COS(i+p)'P-COS(/-p)'P] 


t'mO  p^O  ^ 


3  B 

-lpc22Gn2£-  E  0.,:c  -i^lCOS0-p)'J'-CG50+p)T] 


^  1*0  ^ 


3  ^  ^  B 

-lpc^2Qa^^  x;  0,,^  £^[5//V(i^P)T*5//V(/-p)'P] 


'2  (•>  2 
^  imo  p^o  ^ 


lpc^2QiAU)iGU)R  ^  ^  [S//V(l  +  1)'P-S/^(|-1)'F] 


0; 


-lpc22G(>it/)(Gt^)/?  E  ^  [C0S(i-l)'P-C(?S(/  +  l)4' 


-Ipc^  2g  (AU)R^  ^  ^  Q_^  [C0S(/+p-l)4'  +  C05(i-p  +  1)'F] 


10 


-Ipc^  26  (AU)R^  ^  0.^52  ^hf.lCOSii-^p*\)'V  *  C0S(z-|)-1)'F] 


imo  pmo 


-Ipc^  26  iAU)R^  ^  PPpj  [s/A?(i+p  +  l)4'  -  S//V(i-p-1)'F] 

2  2  ;_y>  »!•/»  4 


l«0  pmo 


Ipc^  26  iAU)R^  ^  ^  [5//V(i+p-l)»p  -  5/Ar(j-p  +  l)'P] 


tmO  p^O 


-Ipc^  26  {AU)R'^  ^  [s/^(i+p-l)4'  +  S/^Cz-p  +  l)*?] 


tmo  p*0 


‘2  fC  W 

-Lpc'^  26  iAU)R  Q  Q_^  j>  [5//y^(,-+p  +  l)4/  +  SlN{i-p-\)^>] 

2  2  p^g 

2  ^  ^ 

-Ipc^  IRS^El^  n  Y,  ®4^ E  [cos(i-p-i)4'  -  cos(i+p  +  i)'p] 

2  2 

-Ipc^  26  (AU)R^  Q  [COSO-p  +  D'F  -  COS(/+p-l)'P].  (29) 

2  2  4 


In  addition,  there  is  an  effect  in  torsion  when  the  pitch  and  elastic  axes  are  different  than  quarter 
chord,  where  the  aerodynamic  force  acts.  The  torsional  moment  due  to  the  pitch  and  elastic  offset 


axes  IS  given  as 


dr  4  ^  dr  dr  Ut 


'£-PA^  ea]^pc  (a  +  Co)  (ef/^  -  CrUp) 
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Performing  the  radial  integration,  again  separating  hannonics  and  letting  (04  -  PA  ej  =  CB, 


=  -CB^pc(fl  +  Co)  +  £0,,  cosin' 


-CB|pc(fl^Co)  £0,,  Sitin' 


-CB^  pc(a  +  Co)i^^  MC)  5^  0.^  [sin  (/♦un'  -  s/yv  (i-1)  *p] 


-  CbI  pc  (a  +  Cfl)  ^■—  (>40)  5^  0„  [cos  (i-l)4'  -  COS  (i  +  1)  S'] 

2  2  imy% 


+  CbI  pc  (AU)^R  5^  0,^  [cos  0  +2)4'  +  COS  (i-2)  4'] 

2  ^  /■a 


^1  (^•'•^ol  1 

+  CbI  pc  2 _ ^ 


1  pc  i _ El  (AU)^R  T  0,,  [sin  0  +  2)4'  +  S//V  (i-2)  4'] 

2  4 


+  CB^  pc  (a*Cp)  ^  (GO) 

_2  ^ 

+  CBlpc (a  +  Co)B3iL  x:  ‘(-Pic  s/iVi4'  +  P,., Cosin') 

^  ^  ino 

+  CbI  pc  {a*Cp)  (AU)  (GO)  B  SW  4' 


CbI  pc  (fl  +  C^)  (AU)R'^  —  Y,  [cos  (i-l)4'  -  COS  (i  +  un*] 


2  •  2 

■^1-0  ^ 


K  .-Q 

CbI  pc  (a  +  Cp)  (AU)  R^  —  Y  —  (i  +  l)4'  -  SIN  (i-l)4']  . 


2-2 
+  1*0  + 
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The  rotor  blade  fli4>ping  equation  of  motion  is 


/a  ^  +  Cb  COS  I'P  +  //,,  SIN  ,'F  .  (33) 

dt  /■o 

where  tg  -  2nd  moment,  Cg  =  flap  damping.  Kg  -  flap  spring,  and  and  -  flap  moment 
forcing  (cosine  and  sine). 


♦  ‘^Pi*  • 

/■o 


Substituting  in  Equation  3  and  equating  5/A^  and  COS  tenns  on  both  sides  of  the  equal  sign. 


and 

and 

Let 

and 


-  /p P,., ^  Cp Pi,  ia*K^ Pi,  -  //i,  (37) 

-  +  <^pP,c  ^pPi,  *  "i,  •  (38) 

.2  ft  ^ 

±P  -  52  i-iCif  Pi,  COS  i'V  -  (-in)2  Pi,  sin  iT  .  (36) 

dt  •■o 

"i.  •  "io-  ‘  Pi.  i-<i.  C39) 

Hi,  -  Hi-  •  Pi,  i-<i,  •  W 


where  =  aero  coefficient  of  the  COS/'f'  terms  and  Pi,(M)i,  =  aero  coefficient  of  the  SIN^f 

terms. 
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Taking  these  terms  over  to  the  left-hand  side  in  Equations  37  and  38  results  in 


-  ^  (Cpifl  -  LAi,)  p,,  +  ATpP,,  -  (41) 

and 

-  /u p„  *  ( -  Cd iQ  -  M„)  ♦  Kp p,,  .  «...  .  (42, 


From  Equations  41  and  42.  the  blade  flapping  responses  p^^  and  p^  can  be  determined. 

Notice  that  the  dynamic  blade  flap  response  feeds  back  into  Equation  6.  This  means  that  the 
blade  flap  response  changes  the  aerodynamic  forcing  function  which  in  turn  means  that  the  dynamic 
blade  flap  motion  will  change.  Therefore,  there  is  nonlinear  coupling  continually  going  on,  and  the 
analysis  should  be  iterated  until  a  solution  is  obtained.  The  analysis  does  indeed  perform  this  iterative 
process  during  the  trim  routine. 

The  dynamic  blade  equations  of  motion  in  torsion  have  the  same  form  as  that  for  flap,  and  elastic 
blade  pitch  occurs,  having  harmonic  terms  and  9^^.  These  elastic  blade  pitch  motions  feed  back 
into  Equation  7  and  also  after  the  aerodynamic  forcing  functions.  In  the  case  of  torsion,  it  is  more 
powerful  than  the  flap  contribution  because  it  gets  directly  into  the  angle  of  attack,  so  that  small 
torsional  deflections  can  cause  large  changes  in  the  blade  aerodynamic  forces. 

4.  APPLICATION  OF  THE  ANALYSIS  FOR  A  UH-60A  HELICOPTER  WITH  A  DAMAGED 
MAIN  ROTOR  BLADE 

For  the  purposes  of  this  report,  the  harmonic  forcing  functions  derived  in  Section  3  are  applied  at 
the  hub  of  the  main  rotor  of  a  UH-60A  helicopter,  and  the  helicopter’s  dynamic  response  is  calculated. 
Vibrations  are  calculated  in  the  cockpit,  assuming  the  fuselage  to  be  a  rigid  body  with  six  degrees  of 
freedom  (i.e.,  longitudinal,  lateral,  vertical,  roll,  pitch,  and  yaw  motions  of  the  helicopter  about  its 
center  of  gravity). 

The  sources  of  the  rotor  blade  and  fuselage  characteristics  data  come  from  a  combination  of 
special  dynamics  data  requested  from  the  manufacturer  and  from  general  technical  manual  data  from 
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the  manufacturer.  There  is  no  one  place  from  which  an  analyst  can  obtain  the  data,  but  it  is  a  matter 
of  digging  and  searching  into  all  kinds  of  documents  or  phoning  the  manufacturer  for  information. 

The  damage  that  is  assumed  is  the  removal  of  sections  of  the  blade  tip  ranging  from  0%  to  48% 
of  the  blade.  The  loss  of  blade  sections  results  in  shonening,  reduction  of  weight,  reduction  of  the 
blade’s  first  and  second  moments,  and  the  increase  of  the  blade’s  first  flap  frequency.  Table  1  shows 
the  changes  of  blade  properties  as  a  function  of  blade  loss.  The  damage  is  applied  to  one  of  the  four- 
bladed  rotor  sets  with  the  helicopter  flying  at  160  knots. 

The  cockpit  vibrations  that  result  from  the  blade  damage  defined  from  Table  1  are  shown  in 
Figure  2.  These  vibrations  are  resultant  vibrations,  the  vector  sum  of  vertical,  longitudinal,  and  lateral 
vibrations. 

It  is  seen  in  Figure  2  that  the  cockpit  vibrations  vary  up  to  2  Gs  at  a  frequency  equal  to  the  rotor 
RPM,  called  1/REV,  in  the  cockpit  when  48%  of  one  blade  is  lost.  Reduced  amounts  of  2/REV, 
3/REV,  4/REV,  and  5/REV  vibrations  are  present  1/REV  means  a  vibration  equal  to  the  rotor  speed, 
2/REV  means  a  vibration  at  two  times  the  rotor  speed,  and  so  on  up  the  line  to  the  higher  /REV 
values. 

5.  DISCUSSION  OF  RESULTS 

These  results  show  the  significance  of  an  unbalanced  rotor  on  vibrations  in  the  cockpit  of  a 
helicopter.  When  a  rotor  is  balanced,  inherent  vibrations  are  transmitted  into  the  fixed  system  fuselage 
at  frequencies  that  are  integer  multiples  of  the  number  of  blades  in  the  rotor.  This  means  for  the  case 
of  the  Black  Hawk  helicopter  with  four  blades,  only  4/REV,  8/REV,  and  12/REV  vibrations  can  be 
transmitted  into  the  fuselage  for  an  undamaged  rotor.  Transmission  of  all  other  frequencies  is  zero. 

For  the  Black  Hawk,  the  4/REV  frequency  with  no  damage  is  .072  Gs,  which  is  insignificant  on 
the  scale  of  vibrations  shown  in  Figure  2.  The  8/REV  frequency  vibrations  are  even  less. 

For  blade  damage,  however,  loss  of  blade  sections  causes  a  steady  unbalanced  load  in  the  rotating 
system,  which  transfers  into  the  fixed  fuselage  system  as  a  1/REV.  When  48%  of  the  outboard  secticm 
of  one  blade  is  removed,  the  vibrations  are  2  Gs,  a  very  severe  condition  for  human  tolerance.  The 
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Table  1.  Rotor  Blade  Physical  Properties 


higher  harmonics  areless  severe  in  magnitude  but  may  be  just  as  effective  in  debilitating  a  human  if 
human  organs  are  resonant  to  these  higher  frequencies  (Linder  1959). 

With  blade  damage,  the  1/REV,  2/REV,  and  3/REV  vibrations  are  all  higher  than  the  ^bient 
4/REV  vibrations  (no  damage). 

The  computer  program  is  an  ongoing  in-house  development  program  conducted  by  the  author. 
Several  papers  have  been  presented  by  the  author  wherein  results  were  generated  using  this  computer 
program  (Fries  1990a,  1990b,  1990c). 

The  computer  program  runs  by  establishing  a  steady-state  trim  condition  and  the  damage  is 
imposed  on  the  main  rotor.  The  vibrations  that  arise  are  instantaneous  in  nature,  of  just  one  rotor 
revolution  duration  after  the  damage  occurs.  The  program  does  not  analyze  the  helicopter  vibrations 
after  the  damage  for  long  periods  of  time.  To  do  this,  we  know  that  the  damage  will  cause  the 
helicopter  to  go  out  of  trim  and  the  pilot  would  have  to  compensate  with  the  controls  in  some  manner. 

The  program  is  of  first  order  accuracy  and  is  considered  to  be  a  first  order  approximation  to  the 
vibrations  that  exist  for  a  helicopter  with  a  damaged  rotor.  The  aerodynamics  are  of  first  order,  the 
second  order  terms  are  neglected,  and  the  airframe  is  a  rigid  body  not  considering  flexibility  or  natural 
modes.  The  program  is  the  first  attempt  at  understanding  this  vibration  problem  and  ongoing  work  is 
needed  to  enhance  the  program.  The  advantage  of  the  program  is  that  it  is  relatively  small  and  easy  to 
use  to  get  quick  results  as  opposed  to  very  large  analyses  existing  at  universities  and  helicopter 
manufacturers  that  require  years  of  experience  and  many  months  to  exercise  a  specific  helicopter 
model. 

The  program  has  approximately  7(X)  lines  of  code  in  BASIC  language  that  runs  on  the  author’s 
Northgate  Elegance  -  425i  personal  computer.  It  is  available  from  the  author  to  anyone  who  wishes 
to  use  it. 
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6.  CONCLUSIONS 


A  woilcable  vibration  analysis  that  can  quantify  aircraft  vibrations  for  a  variety  of  main  rotor 
damages  is  available.  For  this  report,  only  one  type  of  damage  was  considered,  that  of  the  loss  of 
various  amounts  of  the  outboard  sections  of  one  blade. 

This  analysis  brings  up  the  question  of  human  tolerance  to  vibrations  and  the  effect  on  a  pilot’s 
ability  to  fly  the  aircraft  from  a  physiological  point  of  view. 

Another  question  is  the  control  authority  remaining  to  the  pilot  with  blade  loss,  and  then  there  is 
the  further  question  of  the  aircraft’s  stability  characteristics  after  damage. 

This  computer  program  is  a  useful  tool  for  vulnerability  analysis  and  may  be  used  to  investigate 
the  many  factors  that  could  lead  to  helicopter  attrition  caused  by  combat  damage. 
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APPENDIX: 

DESCREPTION  OF  THE  COMPUTER  PROGRAM  TO  PREDICT  HELICOPTER 
VIBRATIONS  FOR  UNDAMAGED  .\ND  DAMAGED  MAIN  ROTOR  BLADES 
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A-1.  DESCRIPTION  OF  THE  RTVIB20  COMPUTER  PROGRAM 


The  computer  code  first  performs  a  trim  analysis  for  forward  speed  level  flight  by  computing  the 
aerodynamic  forces  in  the  manner  described  in  Section  3.  From  these  forces,  a  balance  is  struck  in  the 
vertical  and  horizontal  directions  so  that  level  flight  is  maintained  at  a  constant  forward  speed.  This 
balance  is  maintained  by  adjusting  the  blade  collective  pitch  and  the  aircraft  pitch  attitude  by  iteration 
until  a  balanced  condition  is  obtained. 

The  physical  properties  input  into  the  program  are  listed  in  Section  A-2  of  this  app)endix,  along 
with  the  computer  code  line  numbers.  The  numerical  values  listed  are  for  the  UH-60A  helicopter  and 
come  from  such  a  wide  variety  of  sources  that  they  are  not  listed.  They  come  by  word  of  mouth 
through  telephone  conversations  with  the  manufacturer,  or  by  special  correspondence,  or  by  digging 
out  data  from  the  manufacturer’s  technical  manuals.  Once  a  trimmed  condition  is  achieved,  the 
program  goes  on  to  calculate  fixed  system  vibrations  at  a  position  on  the  fuselage  as  specified  by  the 
user.  In  this  report,  the  cockpit  location  was  chosen. 

Two  iterative  trim  loops  start  at  lines  521  and  530  of  the  code  and  finish  at  lines  1320  and  2350. 
Within  the  inner  loop,  the  blade  collective  pitch  is  varied,  and  in  the  outer  loop  the  aircraft’s  pitching 
attitude  is  changed  until  a  vertical  and  horizontal  force  balance  is  achieved. 

When  the  helicopter  is  flying  straight  and  level,  the  program  will  perform  either  an  undamaged  or 
damaged  main  rotor  vibration  analysis.  Whichever  path  is  taken  depends  on  the  physical  input  data  on 
lines  3060  to  3090.  On  these  input  lines  there  are  inputs  R,  IB,  SB,  W,  A,  CD,  EB,  WN,  ITH,  CM, 
Al,  and  Q  variables  on  N  separate  lines  to  represent  individual  blades.  That  is,  on  line  3060  we  have 
this  data  for  blade  1,  line  3060  for  blade  2,  and  so  on  for  a  four-bladed  rotor.  The  last  line  is  line 
3090.  For  an  undamaged  blade,  the  input  values  on  these  lines  will  be  identical  to  the  corresponding 
values  input  to  lines  50,  300,  360,  440,  or  500,  as  described  in  Section  A-2. 

For  an  undamaged  blade,  all  the  input  parameters  for  one  variable  on  each  of  these  lines 
(3060-3090  for  a  four-bladed  rotor)  are  identical.  When  one  of  the  blades  is  damaged  in  some  way, 
the  physical  properties  will  change  and  this  will  be  reflected  in  a  change  in  one  or  more  of  the  input 
parameters  on  one  of  these  particular  lines.  For  instance,  a  hole  through  one  blade  will  require  IB, 
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SB.  W,  A,  CD,  WN,  ITH,  CM,  Al,  and  Q  on  one  particular  line  for  the  blade  to  be  changed  to 
account  for  damage  (see  input  variable  definitions). 

The  fixed  system  cockpit  vibrations  produced  by  blade  shears  and  moments  for  undamaged  and 
damaged  blades  are  calculated  using  the  computer  code  and  are  reported  in  Section  A-4  for  the 
example  of  the  UH-60A  helicopter. 


A-2.  INPUT  VARIABLES  OF  THE  COMPUTER  CODES 


Variable 

Variable 

Unit  of 

Line 

Numerical 

Incuts 

Inout  Name 

Measurement 

No. 

Inout 

AT 

Pressure  altitude 

ft 

50 

4,000 

T 

Outside  temperature 

deg  F 

50 

95 

CH 

Blade  chord  =  C 

inches 

300 

20.88 

R 

Blade  tip  radius 

ft 

300 

26.83 

TT 

Blade  geometric  twist 

deg/ft 

300 

.671 

OM 

Rotor  speed 

rpm 

300 

258. 

V 

Aircraft  forward  speed 

knots 

300 

160. 

VC 

Vertical  climb  rate 

knots 

300 

0. 

GA 

Rotor  shaft  tilt  =  G, 

+  into  wind 

deg 

300 

3 

PI 

Aircraft  pitch, 

+  into  wind.  Note: 
not  3.14 

deg 

300 

10.49 

CD 

Blade  drag  coefficient, 
steady 

nd 

300 

.04 

GW 

Aircraft  gross  weight 

lb 

360 

16.260. 

TO 

Blade  collective  pitch 

deg 

360 

37.29 

T1 

Blade  lateral  cyclic  pitch 

deg 

360 

.88 

T2 

Blade  longitudinal  cyclic 
pitch 

deg 

360 

-10.79 

N 

Number  of  blades 
per  rotor 

nd 

360 

4 

FE 

Aircraft  profile  flat 
plate  area 

ft^ 

360 

74.55 

IB 

Blade  2nd  flap  moment 

Ib-sec^-ft 

440 

1,570 

SB 

Blade  1st  flap  moment 

lb-sec^ 

440 

87.90 

W 

Blade  weight 

lb 

440 

211 

EB 

Blade  flap  offset  from 
center  of  rotation 

in 

440 

21.98 

PA 

Blade  pitch  axis  from 
leading  edge 

in 

500 

5.22 
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Variable 

Variable 

Unit  of 

Line 

Numerical 

Incuts 

Input  Name 

Measurement 

No. 

Input 

WN 

Blade  1st  torsion  natural 
frequency  nonnalized 
by  rotor  speed 

nd 

500 

3.728 

ITH 

Blade  pitch  inertia 

Ib-sec^-ft 

500 

.865 

CM 

Blade  aerodynamic  pitch 
coefficient  steady 

nd 

500 

-.04 

MC 

Blade  chordwise 
mass  center  from  l.e. 

in 

500 

0 

A1 

Blade  linear  aerodynamic 
pitch  coefficient,  multiplies 
angle  of  attack 

nd 

500 

-.08 

EA 

Blade  elastic  axis 
from  l.e. 

in 

500 

5.22 

Q 

Blade  non-linear 
aerodynamic  pitch 
coefficient,  multiplies 
angle  of  attack  squared 

nd 

500 

0 

lY 

Aircraft  pitch  inertia 

Ib-sec^-ft 

5680 

60,731 

IX 

Aircraft  roll  inertia 

Ib-sec^-ft 

5680 

4,118 

XX 

Longitudinal  distance  from 
rotor  hub  to  the  place 
where  vibrations 
are  calculated 

in 

5680 

134 

YY 

Lateral  distance  to 
vibration  location 

in 

5680 

20 

ZZ 

Vertical  distance  to 
vibration  location 
as  above 

in 

5680 

-100 

EX 

Longitudinal  distance 
from  the  rotor  hub  to 
the  aircraft  c.g. 

in 

5680 

11 

EZ 

Vertical  distance  from 
rotor  hub  to  the  aircraft 

in 

5680 

63 

c.g. 

A-3.  DESCRIPTION  OF  COMPUTER  PROGRAM  OUTPUT 


Only  pertinent  outputs  will  be  given.  There  is  much  output  which  is  intermediate  in  nature  so 
excerpts  will  be  presented. 

After  the  program  iterates  on  a  trim  condition,  a  final  set  of  values  of  blade  collective  pitch, 
aircraft  pitch  attitude,  vertical  force,  and  downwash  velocity  are  produced.  The  results  are  given  in 
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Section  A-4  for  the  undamage  case  for  the  UH-60A  flying  at  160  knots.  Also  presented  in 
Section  A-4  are  the  results  of  calculations  of  the  vertical  fixed  system  hub  shears  and  the  in-plane  hub 
shears.  The  example  case  is  for  an  undamaged  rotor  condition,  in  which  only  the  four/rev  shears  are 
significant  for  a  four-bladed  rotor.  For  a  damaged  case  rotor  blade,  all  the  harmonics  may  be 
significant.  The  1/REV  in-plane  shears  only  are  given,  and  we  see  that  for  the  undamaged  case  these 
shears  are  small;  however,  for  a  damaged  case,  the  steady  damage  in  the  rotating  system  translates  into 
1/REV  fixed  system  shears  that  are  significant 

Fixed  system  4/REV  hub  moments  are  given  in  Section  A-4.  Again,  only  the  4/REV  values  are 
listed  because  only  these  are  transmitted  into  the  fixed  system  for  an  undamaged  rotor.  Also,  in 
Section  A-4  are  listed  the  quasi-theoretical  moments  based  upon  the  expected  output  at  specific 
azimuth  positions.  These  are  used  as  a  check  against  the  values  that  are  calculated  in  the  stepwise 
incrementing  around  the  azimuth.  These  theoretical  moments  come  from  the  harmonic  closed  form 
hand  calculations  from  which  the  magnitudes  at  specific  azimuth  positions  are  known.  For  instance,  at 
1/REV,  we  know  the  cosine  amplitude  at  zero  azimuth  is  either  maximum  or  minimum  and  the  sine 
amplitudes  are  either  maximum  or  minimum  at  the  azimuth  position  of  90  degrees. 

The  final  output  consists  of  cockpit  vibrations  which  are  also  listed  in  Section  A-4.  These  are  the 
base  case  vibrations  for  an  undamaged  rotor  against  which  damage  cases  can  be  compared.  They  are 
broken  down  into  cosine  and  sine  components  of  vibration  in  the  longitudinal,  lateral,  and  vertical 
directions.  Again,  here  only  the  4/REV  vibrations  are  transmitted  from  the  blade  rotating  system  into 
the  fuselage  fixed  system  because  of  a  balanced  rotor  when  no  damage  is  present. 

Note  that  some  of  the  variable  output  names  in  Section  A-4  are  identical  in  physical  meaning,  but 
only  in  a  programming  sense.  That  is,  they  are  variables  in  the  code  that  arc  printed  out  sequentially 
in  time  and  are  reused  at  a  later  point  in  the  execution  of  the  program.  The  line  number  on  which  the 
variable  rests  determines  the  physical  quantity  that  is  printed. 
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A-4.  OUTPUT  VARIABLES  OF  THE  COMPUTER  CODE 


Variable 

Variable 

Unit  of 

Line 

Numerical 

OutDUt 

OutDUt  Name 

Measurement 

No. 

OutDUtS 

TO 

Blade  collective 

deg 

1171 

38.57 

pitch 

1241 

38.57 

PI 

Aircraft  pitch 

deg 

1171 

8.32 

into  wind 

1241 

8.32 

FV* 

Vertical  force 

lb 

1271 

15,948.81 

GW 

Gross  weight 

Ib 

1271 

16,260. 

U 

Downwash  velocity 

ft/sec 

1271 

49.71 

Fixed  System  Vertical  Forces 


I 

(4/REV) 

Harmonic  number 

nd 

4680 

1  to  6 

ACX 

Max.  cosine  amplitude 

lb 

4680 

212.70 

ZCX 

Azimuthal  position 
of  ACX 

deg 

4680 

45 

ACN 

Min.  cosine  amplitude 

lb 

4680 

-212.70 

ZCN 

Azimuthal  position 
of  ACN 

deg 

4680 

0 

ASX 

Max.  sine  amplitude 

lb 

4690 

5.87 

ZSX 

Azimuthal  position 
of  ASX 

deg 

4690 

22.5 

ASN 

Min.  sine  amplitude 

lb 

4690 

-5.87 

ZSN 

Azimuthal  position 
of  ASN 

deg 

4690 

67.5 

In-Plane  Fixed  System  Forces,  1/REV 


FXC 

Last  longitudinal 
cosine  force  calculated 
about  the  azimuthal 
steps 

lb 

4990 

.5019 

ACX 

Max.  longitudinal 
cosine  force 

lb 

4990 

.753 

XCX 

Azimuthal  position  of 
deg  ACX 

deg 

4990 

49.09 

ACN 

Min.  longitudinal 
cosine  force 

lb 

4990 

-.753 

XCN 

Azimuthal  position 
of  ACN 

deg 

4990 

212.7 

♦FV  = 

N*TH*COS(GA-i-PI) 
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Variable 

Variable 

Unit  of 

Line 

Numerical 

Output 

Output  Name 

Measurement 

No. 

Outputs 

FYS 

Last  lateral  sine 

lb 

5000 

-.585 

ASX 

force  calculated 
about  the  azimuthal 
steps 

Max.  lateral  sine 

lb 

5000 

.753 

XSX 

force 

Azimuthal  position 

deg 

5000 

130.9 

ASN 

of  ASX 

Min.  lateral  sine 

lb 

5000 

-.765 

XSN 

force 

Azimuthal  position 

deg 

5000 

310.9 

I 

of  ASN 

Harmonic  number 

Fixed  System  Moments 

nd 

5520 

1  to  6 

(4/REV) 

XSX 

Max.  sine 

ft-lb 

5520 

1.014.56 

PSX 

rolling  moment 
Azimuth  position 

deg 

5520 

67.5 

XSN 

of  XSX 

Min.  sine 

ft-lb 

5520 

-1,014.56 

PSN 

rolling  moment 
Azimuth  position 

deg 

5520 

22.5 

MXS 

of  XSN 

Last  calculated 

ft-lb 

5520 

-.0081 

XCX 

sine  moment  about 
the  azimuth 

Max.  cosine 

ft-lb 

5530 

2.485.99 

PCX 

rolling  moment 
Azimuth  position 

deg 

5530 

0 

XCN 

of  XCX 

Min.  cosine 

ft-lb 

5530 

-2.485.99 

PCN 

rolling  moment 
Azimuth  position 

deg 

5530 

45 

MXC 

of  XCN 

Last  calculated 

ft-lb 

5530 

2.485.99 

YSX 

cosine  moment  about 
the  azimuth 

Max.  sine 

ft-lb 

5540 

2.611.93 

QSX 

pitching  moment 
Azimuth  position 

deg 

5540 

22.5 

YSN 

of  YSX 

Min.  sine 

ft-lb 

5540 

-2,611.93 

pitching  moment 
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Variable 

Variable 

Unit  of 

Line 

Numerical 

OutDUt 

OutDUt  Name 

Measurement 

No. 

Outputs 

QSN 

Azimuth  position 

deg 

5540 

67.5 

MYS 

of  YSN 

Last  calculated 

ft-lb 

5540 

.208 

YCX 

sine  pitching 
moment  about  the 
azimuth 

Max.  cosine 

ft-lb 

5550 

1,083.36 

QCX 

pitching  moment 
Azimuth  position 

deg 

5550 

0 

YCN 

of  YCX 

Min.  cosine 

ft-lb 

5550 

-1,083.36 

QCN 

pitching  moment 
Azimuth  position 

deg 

5550 

45 

MYC 

of  YCN 

Last  calculated 

ft-lb 

5550 

1,083.36 

I 

cosine  pitching 
moment  about 
azimuth 

Theoretical  Fixed  System  Moment  Values 

Harmonic  number  nd  5630 

1  to  9 

(4/REV) 

XMS 

Sine  rolling 

ft-lb 

5630 

-1,014.56 

YMS 

moment 

Sine  pitching 

ft-lb 

5630 

2,611.93 

XMC 

moment 

Cosine  rolling 

ft-lb 

5630 

2,485.99 

YMC 

moment 

Cosine  pitching 

ft-lb 

5630 

1,083.36 

I 

moment 

Harmonic  number 

Cockpit  Vibrations 

nd 

5800 

1  to  10 

(4/REV) 

X2C 

Longitudinal  cosine 

Gs 

5800 

.00140 

X2S 

Longitudinal  sine 

Gs 

5800 

.00412 

Y2C 

Lateral  cosine 

Gs 

5800 

.0578 

Y2S 

Lateral  sine 

Gs 

5800 

-.0235 

Z2C 

Vertical  cosine 

Gs 

5800 

.0367 

Z2S 

Vertical  sine 

Gs 

5800 

.00342 
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Intentionally  left  blank. 


No.  of 

Copies  Organization 


No.  of 

Copies  Organization 


2  Administrator 

Defense  Technical  Info  Center 
ATTN:  DTIC-DDA 
Cameron  Station 
Alexandria,  VA  22304-6145 

1  Commander 

U.S.  Army  Materiel  Command 
ATTN:  AMCAM 
5001  Eisenhower  Ave. 

Alexandria,  VA  22333-0001 

1  Director 

U.S.  Army  Research  Laboratory 
ATTN:  AMSRL-D 
2800  Powder  Mill  Rd. 

Adelphi,  MD  20783-1145 

1  Director 

U.S.  Army  Research  Laboratory 
ATTN:  AMSRL-OP-CI-AD, 

Tech  Publishing 
2800  Powder  MUl  Rd. 

Adelphi.  MD  20783-1145 

2  Commander 

U.S.  Army  Armament  Research. 

Development,  and  Engineering  Center 
ATTN:  SMCAR-IMI-I 
Picatinny  Arsenal,  NJ  07806-5000 


1  Commander 

U.S.  Army  Missile  Command 
ATTN:  AMSMI-RD-CS-R  (DOC) 

Redstone  Arsenal,  AL  35898-5010 

1  Commander 

U.S.  Army  Tank-Automotive  Command 
ATTN:  ASQNC-TAC-DIT  (Technical 
Information  Center) 

Warren.  MI  48397-5000 

1  Director 

U.S.  Army  TRADOC  Analysis  Command 
ATTN:  ATRC-WSR 

White  Sands  Missile  Range,  NM  88002-5502 

1  Commandant 

U.S.  Army  Field  Artillery  School 

ATTN:  ATSF-CSI 

Ft.  Sill,  OK  73503-5000 

(ciMi.  only)  1  Commandant 

U.S.  Army  Infantry  School 
ATTN:  ATSH-CD  (Security  Mgr.) 

Fort  Benning,  GA  31905-5660 

(Unctus.  only) }  Commandant 

U.S.  Army  Infantry  School 
ATTN:  ATSH-CD-CSO-OR 
Fort  Benning,  GA  31905-5660 


2  Commander 

U.S.  Army  Armament  Research, 

Development,  and  Engineering  Center 
ATTN:  SMCAR-TDC 
Picatinny  Arsenal,  NJ  07806-5000 

1  Director 

Benct  Weapons  Laboratory 
U.S.  Army  Armament  Research, 
Development,  and  Engineering  Center 
ATTN:  SMCAR-CCB-TL 
Watervliet,  NY  12189-4050 

(lincbss.  only)  1  Commander 

U.S.  Army  Rock  Island  Arsenal 
ATTN:  SMCRI-IMC-RTTechnical  Library 
Rock  Island,  IL  61299-5000 


WL/MNOI 

EglinAFB,  FL  32542-5000 

Aberdeen  Proving  Ground 

Dir,  USAMSAA 
ATTN:  AMXSY-D 

AMXSY-MP,  H.  Cohen 

Cdr,  USATECOM 
ATTN:  AMSTE-TC 

Dir.  ERDEC 
ATTN:  SCBRD-RT 


Cdr,  CBDA 
ATTN:  AMSCB-CI 


1  Director  1  Uir,  USARL 

U.S.  Army  Aviation  Research  ATTN:  AMSRL-SL-I 

and  Technology  Activity 

ATTN:  SAVRT-R  (Library)  10  Dir,  USARL 

M/S  219-3  ATTN:  AMSRL-OP-CI-B  (Tech  Lib) 

Ames  Research  Center 
Moffett  Field,  CA  94035-1000 
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No.  of 
Copies 

3 


2 


1 

1 

1 

2 


1 

1 


Organization 

Commander 

USA  Aviation  Systems  Command 
ATTN;  AMSAV-ESC,  G.  Kovacs 
SAFE-AV-BH-T,  R.  Olson 
SAFE-AV-AAH-SA,  D.  Roby 
4300  Goodfellow  Blvd. 

St.  Louis,  MO  63120-1798 

Commando' 

AATD 

ATTN:  AMSAT-R-TF,  G.  Hufstetler 
SAVRT-TY-ASV.  H.  Reddick 
USAATCOM 

FLEustis,  VA  23604-5577 

Commander 

ASC 

ATTN;  ASC/XRM,  G.  Bennett 
AFSC  Headquarters 

Wright-Patterson  AFB,  OH  45433-6503 

Commander 

WRDC 

ATTN:  WUnvST.  M.  Untz 
AFSC  Headquarters 

Wright-Patterson  AFB,  OH  45433-6553 
Commander 

ATTN:  ASC/XRYA,  K.  McArdle 
EglinAFB,FL  32542 

Commander 

Naval  Surface  Warfare  Center 
ATTN:  CodeG-13, 

T.  Wasmund 
D.  Dickinson 
Dahlgren,  VA  22448 

McDonnell  Douglas  Helicopter  Company 
ATTN:  William  Sims 
5000  E.  McDowell  Rd. 

Mesa,  AZ  85205 

Boeing  Helicopters  Division 
ATTN:  N.  Caravasos 
P.O.  Box  16858 
Philadelphia,  PA  19142 


No.  of 

Copies  Organization 

1  Silkorsky  Aircraft 
Mail  Stop,  ZIOIA 
ATTN:  G.  Burblis 
6900  Main  Street 
Stratford,  CT  06601-1381 


Aberdeen  Proving  Ground 

1  Dir,  USAMSAA 

ATTN:  AMXSY-AD,  C.  Alston 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laborator)'  undertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes, 
comments/answers  to  the  items/questions  below  will  aid  us  in  our  efforts. 


1.  ARL  Report  Number 


ARL-TR-75 


February  1993 
Date  of  Report _ 


Your 


2.  Date  Report  Received _ 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest  for 

which  the  report  will  be  used.) _ 


4.  Specifically,  how  is  the  report  being  used?  (Information  source,  design  data,  procedure,  source  of 
ideas,  etc.)  _ _ _ 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars  saved, 
operating  costs  avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate. _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate 
changes  to  organization,  technical  content,  format,  etc.) _ 


Organization 


CURRENT  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  Current  or  Correct  address 
above  and  the  Old  or  Incorrect  address  below. 


Organization 


OLD  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 


(Remove  this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  mail.) 


Department  of  the  Army 


OFRCIAL  BUSINESS 


BUSINESS  REPLY  iVLAIL 

FIRST  CUSS  ftaWIT  It  Ml.  HW.  MD 

Postage  *ill  be  paid  by  addiessee 


Director 

U.S.  Army  Research  Laboratory 
ATTN:  AMSRL-OP-CI-B  (Tech  Lib) 
Aberdeen  Proving  Ground.  MD  21005-5066 


